*TROPES

*###
*###

*STUDY 1 (ULINK PHASE I, Survey 2)

*get data
clear
import delimited "/Users/caseyklofstad/Dropbox/u-link phase i survey 2/data/ulink_p1_s2_useable_data.csv"

*select vars to use
keep q1 q5 q4 q3 q2_1 q2_2 q2_3 q2_4 q2_5 q2_6 q11a ///
q22_4 q22_1 q22_2 q22_3 q52 q53 q54 q55 q56 q57 q58 ///
q109 q110 q111 q112 q113 q114 q115 ///
q116 q117 q118 q119 q120 q121 q122 ///
q123 q124 q125 q126 q127 q128 q129 ///
q130 q131 q132 q133 q134 q135 q136 ///
q137 q138 q139 q140 q141 q142 q143

*clean up sample demogs
	*sex
tab q1
generate female=.
replace female=0 if q1==1
replace female=1 if q1==2
tab q1 female
	*age
tab q5
		*convert yob to age
generate age=2020-q5
tab age
			*get median (50th percentile)
sum age, d
	*income
tab q4
	*education
tab q3
	*race
tab q2_1
tab q2_2
tab q2_3
tab q2_4
tab q2_5
tab q2_6
generate nonwhite=.
replace nonwhite=1 if q2_1==.
replace nonwhite=0 if q2_1==1
tab nonwhite
tab q2_1 nonwhite
tab q2_2 nonwhite
tab q2_3 nonwhite
tab q2_4 nonwhite
tab q2_5 nonwhite
tab q2_6 nonwhite

*clean up partisanship
	*1 = dem, 0 = all others
tab q11a
generate dem=0
replace dem=1 if q11a==1
tab q11a dem
label define dem_label 1 "DEM" 0 "all others"
label values dem dem_label
tab q11a dem
	*1 = rep, 0 = all others
tab q11a
generate rep=0	
replace rep=1 if q11a==2
tab q11a rep
label define rep_label 1 "REP" 0 "all others
label values rep rep_label
tab q11a rep

*create the CT scale
	*recode all variables to run low-high on consp beliefs
		*q22_4: Much of our lives are being controlled by plots hatched in secret places
generate plots=.
replace plots=1 if q22_4==5
replace plots=2 if q22_4==4
replace plots=3 if q22_4==3
replace plots=4 if q22_4==2
replace plots=5 if q22_4==1
tab plots q22_4
		*q22_1: Even though we live in a democracy, a few people will always run things anyways
generate few_people=.
replace few_people=1 if q22_1==5
replace few_people=2 if q22_1==4
replace few_people=3 if q22_1==3
replace few_people=4 if q22_1==2
replace few_people=5 if q22_1==1
tab few_people q22_1
		*q22_2: The people who really 'run' the country, are not known to the voters
generate not_known=.
replace not_known=1 if q22_2==5
replace not_known=2 if q22_2==4
replace not_known=3 if q22_2==3
replace not_known=4 if q22_2==2
replace not_known=5 if q22_2==1
tab not_known q22_2
		*q22_3: "big events"
generate big_events=.
replace big_events=1 if q22_3==5
replace big_events=2 if q22_3==4
replace big_events=3 if q22_3==3
replace big_events=4 if q22_3==2
replace big_events=5 if q22_3==1
tab big_events q22_3
	*run PCA analysis
		*note: only 1 factor extracted, explains 68% of variance
factor plots few_people not_known big_events, pcf
rotate
	*get and save PC1 scores as a new variable (factor1)
predict factor1
summarize factor1
	*rescale factor1 to run 0-1 into new variable (consp_beliefs), where 0 = least and 1 = most conspiratorial
generate consp_beliefs=(factor1-(-2.255168))/(1.865602 - (-2.255168))
summarize consp_beliefs
ci means consp_beliefs
		*note: correaltion is r = 1
pwcorr factor1 consp_beliefs, sig

*immigrants replacing
	*create factor variable for treatment
generate immig_repl=.
		*control (q52)
replace immig_repl=1 if q52 !=.
tab q52
tab immig_repl
tab immig_repl q52
		*act now (q53)
replace immig_repl=2 if q53 !=.
tab q53
tab immig_repl
tab immig_repl q53				
		*media elite (q54)
replace immig_repl=3 if q54 !=.
tab q54
tab immig_repl
tab immig_repl q54		
		*common sense (q55)
replace immig_repl=4 if q55 !=.
tab q55
tab immig_repl
tab immig_repl q55		
		*powerful plot (q56)
replace immig_repl=5 if q56 !=.
tab q56
tab immig_repl
tab immig_repl q56		
		*all against us (q57)
replace immig_repl=6 if q57 !=.
tab q57
tab immig_repl
tab immig_repl q57
		*danger if not act (q58)
replace immig_repl=7 if q58 !=.
tab q58
tab immig_repl
tab immig_repl q58		
		*summary
label define immig_repl_label 1 "control" 2 "act now" 3 "media elite" 4 "common sense" 5 "powerful plot" 6 "all against us" 7 "danger if not act"
label values immig_repl immig_repl_label
tab immig_repl
summarize immig_repl
	*clean up DVs: combine into one, recode scale, and flip to run disagree (less conspiratorial) - agree (more conspiratorial)
generate immig_repl_dv=.
replace immig_repl_dv=5 if q52==1 | q53==1 | q54==1 | q55==1 | q56==1 | q57==1 | q58==1
replace immig_repl_dv=4 if q52==6 | q53==6 | q54==6 | q55==6 | q56==6 | q57==6 | q58==6 
replace immig_repl_dv=3 if q52==7 | q53==7 | q54==7 | q55==7 | q56==7 | q57==7 | q58==7 
replace immig_repl_dv=2 if q52==8 | q53==8 | q54==8 | q55==8 | q56==8 | q57==8 | q58==8 
replace immig_repl_dv=1 if q52==9 | q53==9 | q54==9 | q55==9 | q56==9 | q57==9 | q58==9 
tab q52
tab q53
tab q54
tab q55
tab q56
tab q57
tab q58
tab immig_repl_dv
	*OLS
		*base model w/ just treatments
reg immig_repl_dv i.immig_repl
esttab, cells(b(star fmt(2)) se(par fmt(2)))
		*full model w/ controls
reg immig_repl_dv i.immig_repl consp_beliefs dem rep female age q4 q3 nonwhite
esttab, cells(b(star fmt(2)) se(par fmt(2)))
	*ologit as robustness check
		*base model w/ just treatments
ologit immig_repl_dv i.immig_repl
esttab, cells(b(star fmt(2)) se(par fmt(2)))
		*full model w/ controls
ologit immig_repl_dv i.immig_repl consp_beliefs dem rep female age q4 q3 nonwhite
esttab, cells(b(star fmt(2)) se(par fmt(2)))

*deep state
		*create factor variable for treatment
generate deep_state=.
		*control (q109)
replace deep_state=1 if q109 !=.
tab q109
tab deep_state
tab deep_state q109
		*act now (q110)
replace deep_state=2 if q110 !=.
tab q110
tab deep_state
tab deep_state q110
		*media elite (q111)
replace deep_state=3 if q111 !=.
tab q111
tab deep_state
tab deep_state q111		
		*common sense (q112)
replace deep_state=4 if q112 !=.
tab q112
tab deep_state
tab deep_state q112		
		*powerful plot (q113)
replace deep_state=5 if q113 !=.
tab q113
tab deep_state
tab deep_state q113		
		*all against us (q114)
replace deep_state=6 if q114 !=.
tab q114
tab deep_state
tab deep_state q114
		*danger if not act (q115)
replace deep_state=7 if q115 !=.
tab q115
tab deep_state
tab deep_state q115		
		*summary
label define deep_state_label 1 "control" 2 "act now" 3 "media elite" 4 "common sense" 5 "powerful plot" 6 "all against us" 7 "danger if not act"
label values deep_state deep_state_label
tab deep_state
summarize deep_state
	*clean up DVs: combine into one, recode scale, and flip to run disagree - agree
generate deep_state_dv=.
replace deep_state_dv=5 if q109==1 | q110==1 | q111==1 | q112==1 | q113==1 | q114==1 | q115==1
replace deep_state_dv=4 if q109==6 | q110==6 | q111==6 | q112==6 | q113==6 | q114==6 | q115==6 
replace deep_state_dv=3 if q109==7 | q110==7 | q111==7 | q112==7 | q113==7 | q114==7 | q115==7 
replace deep_state_dv=2 if q109==8 | q110==8 | q111==8 | q112==8 | q113==8 | q114==8 | q115==8 
replace deep_state_dv=1 if q109==9 | q110==9 | q111==9 | q112==9 | q113==9 | q114==9 | q115==9 
tab q109
tab q110
tab q111
tab q112
tab q113
tab q114
tab q115
tab deep_state_dv
	*reg
		*base model just treatments
reg deep_state_dv i.deep_state
esttab, cells(b(star fmt(2)) se(par fmt(2)))
		*full model w/ controls
reg deep_state_dv i.deep_state consp_beliefs dem rep female age q4 q3 nonwhite
esttab, cells(b(star fmt(2)) se(par fmt(2)))
	*ologit as robustness check
		*base model just treatments
ologit deep_state_dv i.deep_state
esttab, cells(b(star fmt(2)) se(par fmt(2)))
		*full model w/ controls
ologit deep_state_dv i.deep_state consp_beliefs dem rep female age q4 q3 nonwhite
esttab, cells(b(star fmt(2)) se(par fmt(2)))

*jews control
		*create factor variable for treatment
generate jews_control=.
		*control (q116)
replace jews_control=1 if q116 !=.
tab q116
tab jews_control
tab jews_control q116
		*act now (q117)
replace jews_control=2 if q117 !=.
tab q117
tab jews_control
tab jews_control q117
		*media elite (q118)
replace jews_control=3 if q118 !=.
tab q118
tab jews_control
tab jews_control q118		
		*common sense (q119)
replace jews_control=4 if q119 !=.
tab q119
tab jews_control
tab jews_control q119		
		*powerful plot (q120)
replace jews_control=5 if q120 !=.
tab q120
tab jews_control
tab jews_control q120		
		*all against us (q121)
replace jews_control=6 if q121 !=.
tab q121
tab jews_control
tab jews_control q121
		*danger if not act (q122)
replace jews_control=7 if q122 !=.
tab q122
tab jews_control
tab jews_control q122		
		*summary
label define jews_control_label 1 "control" 2 "act now" 3 "media elite" 4 "common sense" 5 "powerful plot" 6 "all against us" 7 "danger if not act"
label values jews_control jews_control_label
tab jews_control
summarize jews_control
	*clean up DVs: combine into one, recode scale, and flip to run disagree - agree
generate jews_control_dv=.
replace jews_control_dv=5 if q116==1 | q117==1 | q118==1 | q119==1 | q120==1 | q121==1 | q122==1
replace jews_control_dv=4 if q116==6 | q117==6 | q118==6 | q119==6 | q120==6 | q121==6 | q122==6 
replace jews_control_dv=3 if q116==7 | q117==7 | q118==7 | q119==7 | q120==7 | q121==7 | q122==7 
replace jews_control_dv=2 if q116==8 | q117==8 | q118==8 | q119==8 | q120==8 | q121==8 | q122==8 
replace jews_control_dv=1 if q116==9 | q117==9 | q118==9 | q119==9 | q120==9 | q121==9 | q122==9 
tab q116
tab q117
tab q118
tab q119
tab q120
tab q121
tab q122
tab jews_control_dv
	*reg
		*base modle just treatments
reg jews_control_dv i.jews_control
esttab, cells(b(star fmt(2)) se(par fmt(2)))
		*full model w/ controls
reg jews_control_dv i.jews_control consp_beliefs dem rep female age q4 q3 nonwhite
esttab, cells(b(star fmt(2)) se(par fmt(2)))
	*ologit as robustness check
		*base modle just treatments
ologit jews_control_dv i.jews_control
esttab, cells(b(star fmt(2)) se(par fmt(2)))
		*full model w/ controls
ologit jews_control_dv i.jews_control consp_beliefs dem rep female age q4 q3 nonwhite
esttab, cells(b(star fmt(2)) se(par fmt(2)))

*middle east terror via mexico
		*create factor variable for treatment
generate terror_MEX=.
		*control (q123)
replace terror_MEX=1 if q123 !=.
tab q123
tab terror_MEX
tab terror_MEX q123
		*act now (q124)
replace terror_MEX=2 if q124 !=.
tab q124
tab terror_MEX
tab terror_MEX q124
		*media elite (q125)
replace terror_MEX=3 if q125 !=.
tab q125
tab terror_MEX
tab terror_MEX q125	
		*common sense (q126)
replace terror_MEX=4 if q126 !=.
tab q126
tab terror_MEX
tab terror_MEX q126		
		*powerful plot (q127)
replace terror_MEX=5 if q127 !=.
tab q127
tab terror_MEX
tab terror_MEX q127		
		*all against us (q128)
replace terror_MEX=6 if q128 !=.
tab q128
tab terror_MEX
tab terror_MEX q128
		*danger if not act (q129)
replace terror_MEX=7 if q129 !=.
tab q129
tab terror_MEX
tab terror_MEX q129		
		*summary
label define terror_MEX_label 1 "control" 2 "act now" 3 "media elite" 4 "common sense" 5 "powerful plot" 6 "all against us" 7 "danger if not act"
label values terror_MEX terror_MEX_label
tab terror_MEX
summarize terror_MEX
	*clean up DVs: combine into one, recode scale, and flip to run disagree - agree
generate terror_MEX_dv=.
replace terror_MEX_dv=5 if q123==1 | q124==1 | q125==1 | q126==1 | q127==1 | q128==1 | q129==1
replace terror_MEX_dv=4 if q123==6 | q124==6 | q125==6 | q126==6 | q127==6 | q128==6 | q129==6 
replace terror_MEX_dv=3 if q123==7 | q124==7 | q125==7 | q126==7 | q127==7 | q128==7 | q129==7 
replace terror_MEX_dv=2 if q123==8 | q124==8 | q125==8 | q126==8 | q127==8 | q128==8 | q129==8 
replace terror_MEX_dv=1 if q123==9 | q124==9 | q125==9 | q126==9 | q127==9 | q128==9 | q129==9 
tab q123
tab q124
tab q125
tab q126
tab q127
tab q128
tab q129
tab terror_MEX_dv
	*reg
		*base model w/ just treatments
reg terror_MEX_dv i.terror_MEX
esttab, cells(b(star fmt(2)) se(par fmt(2)))
		*full model w/ controls
reg terror_MEX_dv i.terror_MEX consp_beliefs dem rep female age q4 q3 nonwhite
esttab, cells(b(star fmt(2)) se(par fmt(2)))
	*ologit for robustness check
		*base model w/ just treatments
ologit terror_MEX_dv i.terror_MEX
esttab, cells(b(star fmt(2)) se(par fmt(2)))
		*full model w/ controls
ologit terror_MEX_dv i.terror_MEX consp_beliefs dem rep female age q4 q3 nonwhite
esttab, cells(b(star fmt(2)) se(par fmt(2)))

*trump = RUS asset
		*create factor variable for treatment
generate trump_RUS=.
		*control (q130)
replace trump_RUS=1 if q130 !=.
tab q130
tab trump_RUS
tab trump_RUS q130
		*act now (q131)
replace trump_RUS=2 if q131 !=.
tab q131
tab trump_RUS
tab trump_RUS q131
		*media elite (q132)
replace trump_RUS=3 if q132 !=.
tab q132
tab trump_RUS
tab trump_RUS q132
		*common sense (q133)
replace trump_RUS=4 if q133 !=.
tab q133
tab trump_RUS
tab trump_RUS q133	
		*powerful plot (q134)
replace trump_RUS=5 if q134 !=.
tab q134
tab trump_RUS
tab trump_RUS q134		
		*all against us (q135)
replace trump_RUS=6 if q135 !=.
tab q135
tab trump_RUS
tab trump_RUS q135
		*danger if not act (q136)
replace trump_RUS=7 if q136 !=.
tab q136
tab trump_RUS
tab trump_RUS q136		
		*summary
label define trump_RUS_label 1 "control" 2 "act now" 3 "media elite" 4 "common sense" 5 "powerful plot" 6 "all against us" 7 "danger if not act"
label values trump_RUS trump_RUS_label
tab trump_RUS
summarize trump_RUS
	*clean up DVs: combine into one, recode scale, and flip to run disagree - agree
generate trump_RUS_dv=.
replace trump_RUS_dv=5 if q130==1 | q131==1 | q132==1 | q133==1 | q134==1 | q135==1 | q136==1
replace trump_RUS_dv=4 if q130==6 | q131==6 | q132==6 | q133==6 | q134==6 | q135==6 | q136==6 
replace trump_RUS_dv=3 if q130==7 | q131==7 | q132==7 | q133==7 | q134==7 | q135==7 | q136==7 
replace trump_RUS_dv=2 if q130==8 | q131==8 | q132==8 | q133==8 | q134==8 | q135==8 | q136==8 
replace trump_RUS_dv=1 if q130==9 | q131==9 | q132==9 | q133==9 | q134==9 | q135==9 | q136==9 
tab q130
tab q131
tab q132
tab q133
tab q134
tab q135
tab q136
tab trump_RUS_dv
	*reg
		*base model just treats
reg trump_RUS_dv i.trump_RUS
esttab, cells(b(star fmt(2)) se(par fmt(2)))
		*full model w/ controls
reg trump_RUS_dv i.trump_RUS consp_beliefs dem rep female age q4 q3 nonwhite
esttab, cells(b(star fmt(2)) se(par fmt(2)))
	*ologit robustness check
		*base model just treats
ologit trump_RUS_dv i.trump_RUS
esttab, cells(b(star fmt(2)) se(par fmt(2)))
		*full model w/ controls
ologit trump_RUS_dv i.trump_RUS consp_beliefs dem rep female age q4 q3 nonwhite
esttab, cells(b(star fmt(2)) se(par fmt(2)))

*conspiracy against trump
		*create factor variable for treatment
generate against_trump=.
		*control (q137)
replace against_trump=1 if q137 !=.
tab q137
tab against_trump
tab against_trump q137
		*act now (q138)
replace against_trump=2 if q138 !=.
tab q138
tab against_trump
tab against_trump q138
		*media elite (q139)
replace against_trump=3 if q139 !=.
tab q139
tab against_trump
tab against_trump q139
		*common sense (q140)
replace against_trump=4 if q140 !=.
tab q140
tab against_trump
tab against_trump q140	
		*powerful plot (q141)
replace against_trump=5 if q141 !=.
tab q141
tab against_trump
tab against_trump q141		
		*all against us (q142)
replace against_trump=6 if q142 !=.
tab q142
tab against_trump
tab against_trump q142
		*danger if not act (q143)
replace against_trump=7 if q143 !=.
tab q143
tab against_trump
tab against_trump q143		
		*summary
label define against_trump_label 1 "control" 2 "act now" 3 "media elite" 4 "common sense" 5 "powerful plot" 6 "all against us" 7 "danger if not act"
label values against_trump against_trump_label
tab against_trump
summarize against_trump
	*clean up DVs: combine into one, recode scale, and flip to run disagree - agree
generate against_trump_dv=.
replace against_trump_dv=5 if q137==1 | q138==1 | q139==1 | q140==1 | q141==1 | q142==1 | q143==1
replace against_trump_dv=4 if q137==6 | q138==6 | q139==6 | q140==6 | q141==6 | q142==6 | q143==6 
replace against_trump_dv=3 if q137==7 | q138==7 | q139==7 | q140==7 | q141==7 | q142==7 | q143==7 
replace against_trump_dv=2 if q137==8 | q138==8 | q139==8 | q140==8 | q114==8 | q142==8 | q143==8 
replace against_trump_dv=1 if q137==9 | q138==9 | q139==9 | q140==9 | q141==9 | q142==9 | q143==9 
tab q137
tab q138
tab q139
tab q140
tab q141
tab q142
tab q143
tab against_trump_dv
	*reg
	*ologit robusness check
		*base model w/ just treats
reg against_trump_dv i.against_trump
esttab, cells(b(star fmt(2)) se(par fmt(2)))
		*full model w/ controls
reg against_trump_dv i.against_trump consp_beliefs dem rep female age q4 q3 nonwhite
esttab, cells(b(star fmt(2)) se(par fmt(2)))	
	*ologit robusness check
		*base model w/ just treats
ologit against_trump_dv i.against_trump
esttab, cells(b(star fmt(2)) se(par fmt(2)))
		*full model w/ controls
ologit against_trump_dv i.against_trump consp_beliefs dem rep female age q4 q3 nonwhite
esttab, cells(b(star fmt(2)) se(par fmt(2)))


*###

*STUDY 2 (U-LINK PHASE II, Survey 1)

*get data
clear
import delimited "/Users/caseyklofstad/Dropbox/u-link phase ii survey 1/data/U-Link+Phase+II,+Survey+1+(fall+2020)_October+23,+2020_08.34.csv"

*select vars
keep q2 q6 q4 q5 q3_1 q3_2 q3_3 q3_4 q3_5 q3_6 ///
q7 q16_4 q16_1 q16_2 q16_3 ///
q30 q31 q32 q33 q34 q35 q36 ///
q37 q38 q39 q40 q41 q42 q43 ///
q44 q45 q46 q47 q48 q49 q50 ///
q51 q52 q53 q54 q55 q56 q57 ///
q58 q59 q60 q61 q62 q63 q64

*clean up sample demogs
	*sex
tab q2
generate female=.
replace female=0 if q2==1
replace female=1 if q2==2
tab q2 female
	*age
tab q6
	*convert yob to age
generate age=2020-q6
tab age
		*get median (50th percentile)
sum age, d
	*income
tab q4
	*education
tab q5
	*race
tab q3_1
tab q3_2
tab q3_3
tab q3_4
tab q3_5
tab q3_6
generate nonwhite=.
replace nonwhite=1 if q3_1==.
replace nonwhite=0 if q3_1==1
tab q3_1 nonwhite
tab q3_2 nonwhite
tab q3_3 nonwhite
tab q3_4 nonwhite
tab q3_5 nonwhite
tab q3_6 nonwhite
tab nonwhite

*clean up partisanship
	*1 = dem, 0 = all others
tab q7
generate dem=0
replace dem=1 if q7==1
tab q7 dem
label define dem_label 1 "DEM" 0 "all others"
label values dem dem_label
tab q7 dem
	*1 = rep, 0 = all others
tab q7
generate rep=0	
replace rep=1 if q7==2
tab q7 rep
label define rep_label 1 "REP" 0 "all others
label values rep rep_label
tab q7 rep

*create the CT scale
	*recode all variables to run low-high on consp beliefs
		*q16_4: Much of our lives are being controlled by plots hatched in secret places
generate plots=.
replace plots=1 if q16_4==5
replace plots=2 if q16_4==4
replace plots=3 if q16_4==3
replace plots=4 if q16_4==2
replace plots=5 if q16_4==1
tab plots q16_4
		*q16_1: Even though we live in a democracy, a few people will always run things anyways
generate few_people=.
replace few_people=1 if q16_1==5
replace few_people=2 if q16_1==4
replace few_people=3 if q16_1==3
replace few_people=4 if q16_1==2
replace few_people=5 if q16_1==1
tab few_people q16_1
		*q16_2: The people who really 'run' the country, are not known to the voters
generate not_known=.
replace not_known=1 if q16_2==5
replace not_known=2 if q16_2==4
replace not_known=3 if q16_2==3
replace not_known=4 if q16_2==2
replace not_known=5 if q16_2==1
tab not_known q16_2
		*q16_3: "big events"
generate big_events=.
replace big_events=1 if q16_3==5
replace big_events=2 if q16_3==4
replace big_events=3 if q16_3==3
replace big_events=4 if q16_3==2
replace big_events=5 if q16_3==1
tab big_events q16_3
	*run PCA analysis
		*note: only 1 factor extracted, explains 68% of variance
factor plots few_people not_known big_events, pcf
rotate
	*get and save PC1 scores as a new variable (factor1)
predict factor1
summarize factor1
	*rescale factor1 to run 0-1 into new variable (consp_beliefs), where 0 = least and 1 = most conspiratorial
generate consp_beliefs=(factor1-(-2.185178))/(2.006304 - (-2.185178))
summarize consp_beliefs
ci means consp_beliefs
		*note: correaltion is r = 1
pwcorr factor1 consp_beliefs, sig

*chemtrail
	*create factor variable for treatment
generate chemtrail=.
		*control (q30)
replace chemtrail=1 if q30 !=.
tab q30
tab chemtrail
tab chemtrail q30
		*act now (q31)
replace chemtrail=2 if q31 !=.
tab q31
tab chemtrail
tab chemtrail q31				
		*media elite (q32)
replace chemtrail=3 if q32 !=.
tab q32
tab chemtrail
tab chemtrail q32		
		*common sense (q33)
replace chemtrail=4 if q33 !=.
tab q33
tab chemtrail
tab chemtrail q33		
		*powerful plot (q34)
replace chemtrail=5 if q34 !=.
tab q34
tab chemtrail
tab chemtrail q34	
		*all against us (q35)
replace chemtrail=6 if q35 !=.
tab q35
tab chemtrail
tab chemtrail q35
		*danger if not act (q36)
replace chemtrail=7 if q36 !=.
tab q36
tab chemtrail
tab chemtrail q36
		*summary
label define chemtrail_label 1 "control" 2 "act now" 3 "media elite" 4 "common sense" 5 "powerful plot" 6 "all against us" 7 "danger if not act"
label values chemtrail chemtrail_label
tab chemtrail
summarize chemtrail
	*clean up DVs: combine into one, recode scale, and flip to run disagree - agree
generate chemtrail_dv=.
replace chemtrail_dv=5 if q30==1 | q31==1 | q32==1 | q33==1 | q34==1 | q35==1 | q36==1
replace chemtrail_dv=4 if q30==6 | q31==6 | q32==6 | q33==6 | q34==6 | q35==6 | q36==6 
replace chemtrail_dv=3 if q30==7 | q31==7 | q32==7 | q33==7 | q34==7 | q35==7 | q36==7 
replace chemtrail_dv=2 if q30==8 | q31==8 | q32==8 | q33==8 | q34==8 | q35==8 | q36==8 
replace chemtrail_dv=1 if q30==9 | q31==9 | q32==9 | q33==9 | q34==9 | q35==9 | q36==9 
tab q30
tab q31
tab q32
tab q33
tab q34
tab q35
tab q36
tab chemtrail_dv
	*reg
		*base model w/ just treats
reg chemtrail_dv i.chemtrail
esttab, cells(b(star fmt(2)) se(par fmt(2)))
		*full model w/ controls
reg chemtrail_dv i.chemtrail consp_beliefs dem rep female age q4 q5 nonwhite
esttab, cells(b(star fmt(2)) se(par fmt(2)))	
	*ologit for robustnes check
		*base model w/ just treats
ologit chemtrail_dv i.chemtrail
esttab, cells(b(star fmt(2)) se(par fmt(2)))
		*full model w/ controls
ologit chemtrail_dv i.chemtrail consp_beliefs dem rep female age q4 q5 nonwhite
esttab, cells(b(star fmt(2)) se(par fmt(2)))

*bulbs
		*create factor variable for treatment
generate bulbs=.
		*control (q37)
replace bulbs=1 if q37 !=.
tab q37
tab bulbs
tab bulbs q37
		*act now (q38)
replace bulbs=2 if q38 !=.
tab q38
tab bulbs
tab bulbs q38
		*media elite (q39)
replace bulbs=3 if q39 !=.
tab q39
tab bulbs
tab bulbs q39		
		*common sense (q40)
replace bulbs=4 if q40 !=.
tab q40
tab bulbs
tab bulbs q40		
		*powerful plot (q41)
replace bulbs=5 if q41 !=.
tab q41
tab bulbs
tab bulbs q41		
		*all against us (q42)
replace bulbs=6 if q42 !=.
tab q42
tab bulbs
tab bulbs q42
		*danger if not act (q43)
replace bulbs=7 if q43 !=.
tab q43
tab bulbs
tab bulbs q43
		*summary
label define bulbs_label 1 "control" 2 "act now" 3 "media elite" 4 "common sense" 5 "powerful plot" 6 "all against us" 7 "danger if not act"
label values bulbs bulbs_label
tab bulbs
summarize bulbs
	*clean up DVs: combine into one, recode scale, and flip to run disagree - agree
generate bulbs_dv=.
replace bulbs_dv=5 if q37==1 | q38==1 | q39==1 | q40==1 | q41==1 | q42==1 | q43==1
replace bulbs_dv=4 if q37==6 | q38==6 | q39==6 | q40==6 | q41==6 | q42==6 | q43==6 
replace bulbs_dv=3 if q37==7 | q38==7 | q39==7 | q40==7 | q41==7 | q42==7 | q43==7 
replace bulbs_dv=2 if q37==8 | q38==8 | q39==8 | q40==8 | q41==8 | q42==8 | q43==8 
replace bulbs_dv=1 if q37==9 | q38==9 | q39==9 | q40==9 | q41==9 | q42==9 | q43==9 
tab q37
tab q38
tab q39
tab q40
tab q41
tab q42
tab q43
tab bulbs_dv
	*reg
		*base model just treats
reg bulbs_dv i.bulbs
esttab, cells(b(star fmt(2)) se(par fmt(2)))
		*full model w/ controls
reg bulbs_dv i.bulbs consp_beliefs dem rep female age q4 q5 nonwhite
esttab, cells(b(star fmt(2)) se(par fmt(2)))
	*ologit for robustness check
		*base model just treats
ologit bulbs_dv i.bulbs
esttab, cells(b(star fmt(2)) se(par fmt(2)))
		*full model w/ controls
ologit bulbs_dv i.bulbs consp_beliefs dem rep female age q4 q5 nonwhite
esttab, cells(b(star fmt(2)) se(par fmt(2)))
		
*satanic cult
		*create factor variable for treatment
generate satan=.
		*control (q44)
replace satan=1 if q44 !=.
tab q44
tab satan
tab satan q44
		*act now (q45)
replace satan=2 if q45 !=.
tab q45
tab satan
tab satan q45
		*media elite (q46)
replace satan=3 if q46 !=.
tab q46
tab satan
tab satan q46		
		*common sense (q47)
replace satan=4 if q47 !=.
tab q47
tab satan
tab satan q47		
		*powerful plot (q48)
replace satan=5 if q48 !=.
tab q48
tab satan
tab satan q48		
		*all against us (q49)
replace satan=6 if q49 !=.
tab q49
tab satan
tab satan q49
		*danger if not act (q50)
replace satan=7 if q50 !=.
tab q50
tab satan
tab satan q50
		*summary
label define satan_label 1 "control" 2 "act now" 3 "media elite" 4 "common sense" 5 "powerful plot" 6 "all against us" 7 "danger if not act"
label values satan satan_label
tab satan
summarize satan
	*clean up DVs: combine into one, recode scale, and flip to run disagree - agree
generate satan_dv=.
replace satan_dv=5 if q44==1 | q45==1 | q46==1 | q47==1 | q48==1 | q49==1 | q50==1
replace satan_dv=4 if q44==6 | q45==6 | q46==6 | q47==6 | q48==6 | q49==6 | q50==6 
replace satan_dv=3 if q44==7 | q45==7 | q46==7 | q47==7 | q48==7 | q49==7 | q50==7 
replace satan_dv=2 if q44==8 | q45==8 | q46==8 | q47==8 | q48==8 | q49==8 | q50==8 
replace satan_dv=1 if q44==9 | q45==9 | q46==9 | q47==9 | q48==9 | q49==9 | q50==9 
tab q44
tab q45
tab q46
tab q47
tab q48
tab q49
tab q50
tab satan_dv
	*reg
		*base model just treats
reg satan_dv i.satan
esttab, cells(b(star fmt(2)) se(par fmt(2)))
		*full model w/ controls
reg satan_dv i.satan consp_beliefs dem rep female age q4 q5 nonwhite
esttab, cells(b(star fmt(2)) se(par fmt(2)))	
	*ologit for robustness
		*base model just treats
ologit satan_dv i.satan
esttab, cells(b(star fmt(2)) se(par fmt(2)))
		*full model w/ controls
ologit satan_dv i.satan consp_beliefs dem rep female age q4 q5 nonwhite
esttab, cells(b(star fmt(2)) se(par fmt(2)))

*weather control
		*create factor variable for treatment
generate weather=.
		*control (q51)
replace weather=1 if q51 !=.
tab q51
tab weather
tab weather q51
		*act now (q52)
replace weather=2 if q52 !=.
tab q52
tab weather
tab weather q52
		*media elite (q53)
replace weather=3 if q53 !=.
tab q53
tab weather
tab weather q53
		*common sense (q54)
replace weather=4 if q54 !=.
tab q54
tab weather
tab weather q54	
		*powerful plot (q55)
replace weather=5 if q55 !=.
tab q55
tab weather
tab weather q55	
		*all against us (q56)
replace weather=6 if q56 !=.
tab q56
tab weather
tab weather q56
		*danger if not act (q57)
replace weather=7 if q57 !=.
tab q57
tab weather
tab weather q57
		*summary
label define weather_label 1 "control" 2 "act now" 3 "media elite" 4 "common sense" 5 "powerful plot" 6 "all against us" 7 "danger if not act"
label values weather weather_label
tab weather
summarize weather
	*clean up DVs: combine into one, recode scale, and flip to run disagree - agree
		*note: q51 scale 1/6/7/8/9; rest = 1-5
generate weather_dv=.
replace weather_dv=5 if q51==1 | q52==1 | q53==1 | q54==1 | q55==1 | q56==1 | q57==1
replace weather_dv=4 if q51==6 | q52==2 | q53==2 | q54==2 | q55==2 | q56==2 | q57==2 
replace weather_dv=3 if q51==7 | q52==3 | q53==3 | q54==3 | q55==3 | q56==3 | q57==3 
replace weather_dv=2 if q51==8 | q52==4 | q53==4 | q54==4 | q55==4 | q56==4 | q57==4 
replace weather_dv=1 if q51==9 | q52==5 | q53==5 | q54==5 | q55==5 | q56==5 | q57==5 
tab q51
tab q52
tab q53
tab q54
tab q55
tab q56
tab q57
tab weather_dv
	*reg
		*base model just treats
reg weather_dv i.weather
esttab, cells(b(star fmt(2)) se(par fmt(2)))
		*full model w/ controls
reg weather_dv i.weather consp_beliefs dem rep female age q4 q5 nonwhite
esttab, cells(b(star fmt(2)) se(par fmt(2)))
	*ologit for robustness check
		*base model just treats
ologit weather_dv i.weather
esttab, cells(b(star fmt(2)) se(par fmt(2)))
		*full model w/ controls
ologit weather_dv i.weather consp_beliefs dem rep female age q4 q5 nonwhite
esttab, cells(b(star fmt(2)) se(par fmt(2)))

*covid-19 look-alikes
		*create factor variable for treatment
generate alikes=.
		*control (q58)
replace alikes=1 if q58 !=.
tab q58
tab alikes
tab alikes q58
		*act now (q59)
replace alikes=2 if q59 !=.
tab q59
tab alikes
tab alikes q59
		*media elite (q60)
replace alikes=3 if q60 !=.
tab q60
tab alikes
tab alikes q60
		*common sense (q61)
replace alikes=4 if q61 !=.
tab q61
tab alikes
tab alikes q61	
		*powerful plot (q62)
replace alikes=5 if q62 !=.
tab q62
tab alikes
tab alikes q62		
		*all against us (q63)
replace alikes=6 if q63 !=.
tab q63
tab alikes
tab alikes q63
		*danger if not act (q64)
replace alikes=7 if q64 !=.
tab q64
tab alikes
tab alikes q64
		*summary
label define alikes_label 1 "control" 2 "act now" 3 "media elite" 4 "common sense" 5 "powerful plot" 6 "all against us" 7 "danger if not act"
label values alikes alikes_label
tab alikes
summarize alikes
	*clean up DVs: combine into one, recode scale, and flip to run disagree - agree
generate alikes_dv=.
replace alikes_dv=5 if q58==1 | q59==1 | q60==1 | q61==1 | q62==1 | q63==1 | q64==1
replace alikes_dv=4 if q58==2 | q59==2 | q60==2 | q61==2 | q62==2 | q63==2 | q64==2 
replace alikes_dv=3 if q58==3 | q59==3 | q60==3 | q61==3 | q62==3 | q63==3 | q64==3 
replace alikes_dv=2 if q58==4 | q59==4 | q60==4 | q61==4 | q62==4 | q63==4 | q64==4 
replace alikes_dv=1 if q58==5 | q59==5 | q60==5 | q61==5 | q62==5 | q63==5 | q64==5 
tab q58
tab q59
tab q60
tab q61
tab q62
tab q63
tab q64
tab alikes_dv
	*reg
		*base model just treats
reg alikes_dv i.alikes
esttab, cells(b(star fmt(2)) se(par fmt(2)))
		*full model w/ controls
reg alikes_dv i.alikes consp_beliefs dem rep female age q4 q5 nonwhite
esttab, cells(b(star fmt(2)) se(par fmt(2)))
	*ologit for robustness
		*base model just treats
ologit alikes_dv i.alikes
esttab, cells(b(star fmt(2)) se(par fmt(2)))
		*full model w/ controls
ologit alikes_dv i.alikes consp_beliefs dem rep female age q4 q5 nonwhite
esttab, cells(b(star fmt(2)) se(par fmt(2)))
